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Abstract 

In a recent preprint by Deutsch et al. [5] the authors suggest the pos- 
sibility of polynomial approximability of arbitrary unitary operations on n 
qubits by 2-qubit unitary operations. We address that comment by proving 
strong lower bounds on the approximation capabilities of g-qubit unitary 
operations for fixed g. We consider approximation of unitary operations 
on subspaces as well as approximation of states and of density matrices by 
quantum circuits in several natural metrics. The ability of quantum circuits 
to probabilistically solve decision problem and guess checkable functions is 
discussed. We also address exact unitary representation by reducing the 
upper bound by a factor of n 2 and by formalizing the argument given by 
Barenco et al. [1] for the lower bound. The overall conclusion is that almost 
all problems are hard to solve with quantum circuits. 

1 Introduction 

There has recently been great interest in the properties of quantum computa- 
tion and quantum circuits, partly due to the proof by Shor [9] that quantum 
computation can be used to efficiently factor and compute discrete logarithms 
in modular arithmetic. An important problem in quantum computation is to 
better elucidate the relationship between quantum complexity and classical com- 
plexity. It is well known that a quantum computation can be simulated by a 

'Address: MS B265, Los Alamos National Laboratory, Los Alamos, NM 87545. This work 
was performed under the auspices of the U.S. Department of Energy under Contract No. 
W-7405-ENG-36. 

1 



in 
o\ 
o\ 

< 
00 

O 

o 

00 

O 

o> 

Oh 



classical Turing machine with at most an exponential speedup. In addition to 
Shor's result, there is an accumulating body of evidence to support the idea 
that this exponential speedup is necessary. This includes works by Deutsch [6], 
Bernstein and Vazirani [4] and Simon [10], which show that relative to cer- 
tain oracles quantum computation is exponentially more efficient then classical 
computation, even in the probabilistic model. 

Here we address some open issues in quantum circuit complexity. In 
classical circuit theory a simple counting argument shows that almost all func- 
tions from 2^ (the space of ra-bit patterns) to 2^ are hard in the sense that to 
compute them requires exponentially many gates. In quantum circuit theory, a 
similar result is known to hold for exact computation. However, the question 
of an exponential lower bound for approximate or probabilistic computation 
of both classical and non-classical problems had not been answered. In fact, 
in a recent paper by Deutsch, Barenco and Ekert [5], an intuitive argument is 
given for why they believe that approximate computation with quantum circuits 
might be possible with polynomially many two-bit gates. In this paper we show 
that this is not the case. In the process we investigate the overall difficulty of 
various related computation tasks in the quantum circuit model. 

We present a general approach to studying worst case quantum circuit 
complexity by considering the problems of representing and approximating ac- 
tions of unitary operators on arbitrary sets of orthonormal vectors in the Hilbert 
space. We consider approximation of states in the usual norm as well as approx- 
imation of density matrices in a fairly weak norm. For the former problem we 
show that unless the number of gates is exponentially large, almost all states are 
approximated no better than random by any circuit. For the second problem 
we use the total variation distance and obtain a similar (though suboptimal) 
result. 

We define two types of classical problems, one which requires a proba- 
bilistic solution and another which requires a checkable solution. In each case we 
give explicit bounds on the measure or number of problems that can be solved 
by quantum circuits with a given number of gates. The bounds all imply that 
almost all problems are exponentially hard for quantum circuits. 

An interesting consequence of this work is that even the approximate 
state conversion problem is exponentially difficult. That is, for almost all states 
(and almost all pairs of states and density matrices), the smallest quantum 
circuit which converts the first state to an approximation of the second requires 
exponentially many gates. This should be compared to the classical theory of 
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relative information distance [3]. It reflects the additional information that is 
implicitly available in quantum states. 

An outline of the paper follows. Section 2 summarizes the notation 
and quantum computation background required for this work. In Section 3, we 
revisit the upper and lower bounds for exact unitary representation given in [1], 
improving the upper bound by a factor of n 2 and proving the lower bound 
in [1]. Section 4 covers approximation of the action of unitary operators on 
members of a basis in the usual Hilbert space norm. We also obtain bounds for 
approximation of density matrices from the point of view of measurement. In 
Section 5 we consider two types of classical problems and count the number of 
such problems that cannot be solved with bounded size quantum circuits. 

2 Preliminaries 

A basic understanding of Hilbert spaces and measure theory is assumed. The 
k xk identity matrix is denoted by 1\. or simply / if k is understood. A diagonal 
matrix with diagonal elements given by d = ( d\, . . . , dj~ ) is denoted by D(d). 

Quantum circuits. To understand the problems in quantum circuit complex- 
ity requires some knowledge of quantum circuits and their behavior. Descrip- 
tions of quantum circuits can be found in [12] and [1]. The input space of 
quantum circuits operating on n qubits is a complex Hilbert space of dimension 
2 n which is represented as an ra-fold tensor product Q® n of the two dimen- 
sional qubit space Q generated by the orthonormal basis {|0), |1)} (using Dirac 
notation). For example, Q® 2 has the standard orthonormal basis 

{|0)|0),|0)|1),|1)|0),|1)|1)}, 

where juxtaposition of vectors denotes tensor product. The standard orthonor- 
mal basis is naturally identified with bit patterns in 2^, and we write \b) for 
|6o)-"|^n-i) 5 where b{ is the i'th bit of bit pattern (or binary number) b. The 
expressions n and l n represent the bit patterns of n O's and n l's, respectively. 
A number k used in the context of a bit pattern denotes it's binary representa- 
tion (in reverse order). In general, we will use \u), \v), . . . for normalized states 
and u,v, . . . to denote arbitrary, not necessarily normalized vectors in Hilbert 
space. 

The functions of quantum computation are unitary operations on the 
input space. The set of unitary operations on Q® n is denoted by U(Q® n ) or 
U n for short. The output of a quantum circuit is therefore an element of Q® n . 
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Note that unitary operations are invertible, so a quantum computation is always 
reversible. Reversibility is not a strong restriction from the point of view of 
complexity theory, as discussed in [2, 8]. 

The correspondence between quantum circuits and classical (reversible) 
circuits is established by restricting the input of a quantum circuit to the ele- 
ments of the standard basis, and by measuring the output of the circuit. The 
measurement usually results in projecting the output state \u) onto one of the 
basis elements \b). The probability of observing \b) as the projected output is 
given by Prob(6 | \u)) = |(m|6)| 2 . Thus the general quantum circuit exhibits 
probabilistic behavior. 

In many cases, the relevant domain of a quantum circuit is restricted to 
states of the form \bO n ~ l ) and instead of measuring all output bits, only the first 
/ bits may be observed. If the output state is \u), the probability of observing 
that the first / bits are b is Prob(6 | \u)) = J2b' \ {u\bb')\ 2 . 

Quantum gates. The notion of quantum circuit complexity depends on the 
circuit elements or quantum gates that can be used to build a unitary operation. 
The most natural definition requires that each quantum gate act locally on the 
input space. This means that the unitary operation corresponding to the gate 
should involve only a small number of qubits. More formally, let S C [n] = 
{1, . . . , n}. Write Qf 71 and Qf n for the tensor product of the factors of Q® n at 
positions in S and not in S , respectively. Then Q® n is naturally isomorphic to 
Qs n ® Q§ n - We can take a unitary operator in U(Q ( § n ) ~ U\ s \ and have it act 
on Q® n by taking its tensor product with the identity. A g-qubit gate acting 
on Q® n is a member of U g acting on Q® n via a g-qubit factor Q®™ ■ y We 
denote the family of g-qubit gates by G g . Note that for g < g', G g C G g i. 

Quantum circuit complexity, like classical circuit complexity, concerns 
the question of what operations can be represented by compositions of b gates. 
The first question to be answered is of course whether the gates chosen are 
universal. This was solved by DiVincenzo and others [7]: Sufficiently large com- 
positions of members of G2 can represent any unitary operation. Let G g b ' n ^ be 
the set of all unitary operations in U n which can be represented as a compo- 
sition of b members of G g . In general, let y( 6 ' n ) be the set of operators in U n 
representable by compositions of b operators in V. 

Metrics for U n . In order to study the complexity of approximating states, 
unitary operations or other types of functions, we need suitable metrics. There 
are several choices and in general we will make use of the weakest one. We 
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consider linear operations as matrices over the standard basis. For A a linear 
operation on Q® n , its matrix is denned by A^ y = {b\A\b'). Most of the metrics 
to be considered are related to the standard Euclidean norm on a Hilbert space 
H defined by \u\ = \Ju*u. 

The Frobenius norm. The Frobenius norm of A is defined by 



\\a\\ f = Eivp 

V b,b' 

= VtrA*A. 
The two-norm. The two-norm of A is given by 
||A||2 = sup |A|a;)|. 

\x) 

The weak two-norm. This is a weak version of the two-norm which 
is useful for quantum computation: 

||A|| 2 ,fc = max|A|6)|. 

Total variation distance. The output of a quantum circuit induces 
a probability distribution on measurement outcomes. For two probability dis- 
tributions fj, and is, the total variation distance is defined by 

v((i,v) = £>(&)- 

b 

For two states \u) and \v), their total variation distance on the first / bits is 
given by 

vi(\u),\v)) = v(Fmb(- | |w)),Prob(- | \v))) 

= J2 |Prob(6 | |ti))-Prob(6 | \v))\. 

Total variation distance for unitary operators. For unitary op- 
erators U and V, we can define a total variation distance by 

v lik (U,V) = m ax .v,(U\b),V\b)). 

The following lemmas summarize the relationships between the differ- 
ent distance measures. 

Lemma 2.1 The following inequalities hold: 

M\ u )A v )) < v i+v(\ u ), \v)), 

vi(\u),\v)) < 2\\u)-\v)\. 
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Proof. We can consider the v\ seminorms as the L\ distance of the probability 
distributions P u {b) = Prob(6 | \u)) and P v (b) = Prob(6 | \v)) considered as 
vectors over b. Let P x ,i(b) = Prob(6i | \x)) for i = 0,1. Then P x = P Xt o + P x ,ii 
which implies that 

Vl(\u),\v)) = \P U -Pv\l 

= \Pu,0 — Pv,0 + Pu,l — Pv,l\l 
< \Pu,0 — Pvfl\l + \Pu,l — Pv,l\l 
= V{ +1 . 

This gives the first inequality. 

To see the other inequality requires going back to the definitions. For a 
vector x or state \x), let Xf, be the projection onto the space spanned by {\bb')}^, 



J2{bb'\x)\bb'). 



Then Prob(6 | \u)) 



|2 



V{(\u), \v)) = \\ U b\ 2 ~ Kl 

b 



£ \/? l "'- t " l2 l\/? l, " l2+ \/? lt " l2 / 

= 2|| M )-|«)|, 

where we used the Schwarz inequality in the second to last step. 

Lemma 2.2 The following inequalities hold: 

\\A\\ 2 ,k < \\A\\ 2 

< \\ a \\f, 

ll^-lb,*; < ll^-lb.fc+fcs 

vi,k(U,V) < vi + i' tk +k'(U,V), 

vi,k(U,V) < 2\\U-V\\ 2 , k . 
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Proof. To see that ||A||2 < write 

\ Au \ 2 = X] I ^Z A b,b'Ub'\ 2 

b b' 

< E(El^Vl 2 )M 2 

b v 

= \\A\\ 2 F \u\ 2 , 

by the Schwarz inequality. The fact that increasing k in the subscript in- 
creases the norm is true by definition. The remaining inequalities follow from 
Lemma 2.1. ■ 

According to Lemmas 2.1 and 2.2, if it is difficult to approximate states 
or operators in a vi or (semi)metric, then all the other metrics are difficult 
too. Note that vi(\u), \v)) can be viewed as comparing the density matrices 
induced by \u) and \v) in the space spanned by the first / bits by comparing 
only the diagonal of the density matrices. 

We will make use of the following lemma which describes the behavior 
of the 2-norm in conjunction with composition of unitary operators. 

Lemma 2.3 Fori £ {1,2}, letUi andVi be unitary operators. If \\Ui — Vi\\ 2t k < 
S and \\U 2 - V 2 \\2,k < P, then \\U-JJ2 - Vi^lb.fc < <*> + P- 

Proof. 

\\U 1 U 2 -V 1 V 2 \\ 2 , k = \\U t U 2 - U t V 2 - (V t V 2 - UtV^k 

< \\U!(U 2 - V 2 )\\ 2 , k + \\iVt - UtW^k 

< S + P- 

where unitarity was used in the last step. ■ 

Measures for U n . To "count" the fraction of unitary operations that can be 
approximated by 6-gate circuits requires suitable measures on U n . We will con- 
sider two measures, depending on whether we are interested in approximating 
states or approximating probability distributions. Both measures are induced 
by measures on sets of k orthonormal states. For a unitary operator U, let 

MU) = ( u\i)\o <i<k). 

The sphere measure. The unit vectors in a Hilbert space H of 
dimension N are elements of the unit sphere S 2 n in 2N real dimensions. The 
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Euclidean metric induces a measure fj, on S2N f° r which 

KS 2 n) = 2 



T(NY 



Let = ^s 2N ) ^ e ^ e normai i ze( i measure. Let Uk be the set of sequences 
u = ( uo, . . . ,Uk-i ) of orthonormal states in H. One can view an element 
u G Uk as being obtained by choosing \uq) on the unit sphere in Q n , then 
choosing \ui) on the induced unit sphere in the subspace orthogonal to \uo) and 
so on. This induces a measure on Uk by integrating functions "inside out" with 
respect to ^1, that is by integrating over the fc'th element of the sequence of 
states first. This measure is denoted by fj,k- For measurable functions /(u), 

cjti fc (u)/(u) = J dfJ,k-i({ u , . . -,u k -2 )) J dfJ,i(uk-i)f(u), 

where fj,i m the inner integral is over a sphere of dimension 2(N — k + 1). The 
measure fj,k induces a measure on U n via <$>\.. We denote this measure by fj,k as 
well. Note that the induced measure on U n is defined on a u-subalgebra of the 
Borel measurable sets. For our purposes, we will extend fj,k to all subsets of U n 
by defining fj,k(X) = p,k(<f>k(X)), where p,k is the outer measure. 

For us, the most important property of fj,i is the measure of the sphere 
within e of a state. For any vector u, let B U)t be the set of vectors v satisfying 
\ u — v \ < e - 

Lemma 2.4 Let \u\ = 1 and let S be a unit sphere in a subspace of complex 
dimension m. Then for e < y/2 and m > 3, 



2m-l 



V2m- 1(1 - e 2 / 2 )' 
The proof of Lemma 2.4 is in the appendix. 

The density measure. Consider the map Prob; which takes a vector 
u to the vector Prob;(M) defined by Prob;(M)fe = J2b' \ u bb'\ 2 where b is an l- 
bit string. For unit vectors u, Prob(M) is in the simplex A(2') of dimension 
2' defined by J2b Prob(M)6 = 1 and for each b, Prob(M)fe > 0. The Euclidean 
metric induces a measure v on A(iV) for which v(A) = p^y- Let = -^^j 
be the normalized measure. This measure induces a measure on Uk via the 
product measure and the map u ( Prob;(Mo), . . . , Prob;(Mfc_i) ). Note that if 
2 n > k2\ then this map is onto so that the induced measure is non-trivial. This 
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again induces a measure on U n , also denoted by z//^, via fa. This measure is 
extended to arbitrary subsets as we did for 

We need an analog of Lemma 2.4 for z/;^. For any state \u), let B u ^ t 
be the set of states \v) satisfying vi(\u), \v)) < e. 

Lemma 2.5 v hl {B uU ) < ^e) 2 '" 1 . 

The proof of Lemma 2.5 is in the appendix. The result can be strength- 
ened for e > ^, but this is not included in this work. 

3 Exact representation of unitary operators 

DiVincenzo [7] showed that G^'^ = U n for 6 sufficiently large. The open prob- 
lem is to determine 

6(2, n) = min{6 | d, M = U n }. 

Barenco et al.'s [1] analysis shows that 6(2, n) = 0(n 3 4 n ) for DiVincenzo's 
construction. We will show that 6(2, n) = 0(n4 n ). They also give an intuitive 
argument for 6(2, n) > 4 n /9 — l/3ra — 1/9, which we formalize below. 

It is interesting to consider a more detailed study of unitary represen- 
tation by gates. For sequences of orthonormal states u = ( |wo), . . . , ) and 
v = ( |fo), . . ., ) let b(g,n; u; v) be the smallest b such that there exists 
U G Gg'' n ^ with U\ui) = \v{) for each i. If the \u{) are the first i standard basis 
vectors, then we will write 6(</,n;v) for 6(</,n;u;v). Note that if U G Gg'' n \ 
then U~ x G G g (b,n), which implies that 

Ha, n i u ; v ) < Ha, n ; v ) + Ha, n ; u )- 

Let b(g, n; k) = max{6(</, n; u)}, where the maximum is over all possible choices 
of the \ui). 

Let I w be the unitary operation defined by I w \0) = e™|0) and I w \b) = 
|6) for 6^0. Another complexity parameter of interest is the following: 

b(g,n;Q = min{6 | Vwl w e Gf' n ^}. 

I w is a special case of the A n -i(U) gates described in [1], where it is shown that 
6(2, n; i*) = 0(n 2 ). The reason for considering b(g,n;I*) is explained by the 
next theorem. 

Theorem 3.1 b(g, n; k) < k(2b(g, n; 1) + b(g, n; /*))• 



9 



Proof. Let u = ( \uo), . . . , ) be k orthonormal states. The proof is 

in two steps. The first step is to show that the transformation which maps 
\i) to \ui) can be extended to a unitary operator with at most k non-identity 
eigenvalues. The second is to observe that by diagonalizing the unitary opera- 
tion and exploiting 1-transitivity, it can be decomposed into a product of i^'s, 
conjugated by operators instantiating 1-transitivity for the eigenvectors. 
The first step is accomplished by the following lemma. 

Lemma 3.2 Let V be a k X to matrix whose rows are orthonormal. Then there 
exists an (to — k) X to matrix Z with orthonormal rows such that the unitary 
matrix obtained by extending V by Z has at least to — k eigenvalues with value 
1. 



Proof. Let W be an arbitrary orthonormal basis of the subspace orthogonal 
to the row space of V. The desired matrix Z will be of the form UW for U 
unitary. The problem is solved if we can determine U and and an to X (to — k) 
matrix X with column space of dimension (to — k) such that 



V 

UW 



X = X. 



To see that there exists such an X satisfying that VX consists of the first k 
rows of X, it suffices to observe that V — [ifc,0] has a null space of dimension 
at least to — k. For this X, write 



X 



x 1 
x 2 



where X\ is k X (to — k). 

We have the following identities: 



x*x = xix-l + x;x 2 



X{X X 

x*w*wx 



1-m — k i 

x*v*vx, 

X*(L m - V*V)X 

Im-k — X*[X\ 

XnXo. 
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The identity (WX)*(WX) = X 2 X 2 implies that there exists a unitary U such 
that UWX = X 2 , as desired (see Lemma A.l in the appendix). ■ 

Extend the \ui) to a orthonormal basis such that the matrix U whose 
columns are the basis states has at most k eigenvalues different from 1. This 
allows us to write U = V(D ® I 2 n_ k )V*, where D = D(e tw °,. . . , e™*- 1 ). Let V t 
be unitary operations representable by b(g,n; 1) 2-qubit gates which take |0) to 
the i'th column of V (the i'th eigenvector of U). Then 

U = V k I Wk V k * ...Vtl^V*. 



To obtain a general upper bound on b(g,n;k) still requires bounding 
b(g,n; 1). We show that 6(2, n; 1) = 0(n2 n ), thus improving the upper bound 
given in [1] by a factor of n 2 for the case k = 2 n . 

Theorem 3.3 6(2, n; 1) = 0(n2 n ). 

Proof. Suppose we would like to map |0) to \u). Write 

\u)= J2 |&>«M)|0) + a(M)|l)). 

6e2["-!] 

Let A 2 (b) = |a(6,0)| 2 + |a(6, 1)| 2 . Then Y,b A2 ( b ) = 1 - The idea is as follows: 
Using induction, apply an operation to the first n—1 bits, making sure that the 
amplitude of |6) is A(b). Next apply a conditional operation on the last bit to 
distribute the amplitude correctly for each |6). 

More formally, let \v) = ^be^^"- 1 ^ Choose a unitary operator 

U on n — 1 qubits which maps |0) to |f)|0) using at most 6(2, n — 1; 1) gates. 
Next choose conditional operations Ub on n qubits which have the property that 
Ub\b') = 1 6') unless the first n — 1 bits of 6' are 6. In that case define 

^|M) = I ^|6)(a(6,0)|0) + a(6,l)|l)). 

The composition of U with each of the Ub has the desired effect on |0). It remains 
to determine the complexity of Ub- Note that we have made no requirement on 
Ub\bl). This means that without loss of generality, det(Ub) = 1. This implies 
that Ub is one of the conditional operators A n -i(Vb) with V& in SU(2) discussed 
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in [1]. It is shown there that A n -i(Vb) can be constructed with 0(n) gates. This 
gives 

6(2, n; 1) < cn2 n ~ 1 + 6(2, n - 1; 1), 
for some constant c. Thus 6(2, n; 1) = 0(n2 n ). ■ 

Lower bounds on b(2,n;k) can be obtained by using dimensional ar- 
guments. We repeat the argument given in [1], with the relevant details filled 
in to obtain a more general result. 

Theorem 3.4 6(2, n; k) > |(2 n+1 - k) - \n - \ . 

Proof. Consider a quantum circuit implementing operator U with the gates 
§i labeled by the corresponding unitary operators JJ{ G Ui acting on two qubits. 
Two circuits are structurally the same if the only difference is in the associated 
unitary operators. Note that there are finitely many structurally distinct circuits 
involving 6 gates (at most Q) 6 )- The set of sequences ( U\0), . . . , U\k — 1) ) 
achievable by a fixed structure S is the range of a smooth function f$ with 
domain the sequence ( . . . ,Ui, . . . ). The range of f$ is to be compared to Uk 
whose dimension is k(2 n+1 — k). Clearly, the dimension of the domain manifold 
is 6dim([/2)- The dimension of U2 is 16. The task is to determine an upper 
bound on the dimension of the range of f$. The following lemmas reduce the 
initial estimate of 166. 

Lemma 3.5 Any gate other than the first one adds at most 15 dimensions to 
the range of S . 

Proof. We can remove the global phase from JJ{ for i > 2 without changing 
the range. In other words, except for Uo, we can restrict the JJ{ to unitary 
operators with at least one fixed point. ■ 

Lemma 3.6 If gate gi satisfies that one of its inputs is an output of another 
gate, then the effective contribution of JJ{ to the dimension of the range of S is 
at most 12. 

Proof. Without loss of generality assume that the input which comes from 
another gate is the first input of gate g;. For every V G U2, there is a neighbor- 
hood of V that can be parametrized in the form XW, where W is an operator 
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which acts only on the first input qubit and X lives in a submanifold of dimen- 
sion 12. This follows from a fundamental theorem for homogeneous manifolds, 
see [11], page 120. The operator W acts only on one output of the previous gate 
and can therefore be absorbed by that gate. In the complete circuit, one can 
work one's way from the output to the input, restricting JJ{ to one of a finite 
number of local neighborhoods covering U2, restricting f$ to the corresponding 
X's and moving the W 7 s to the previous gate. The result is a decomposition of 
the domain of f$ to a finite union of restrictions to smaller dimensional spaces 
without affecting the union of the ranges. ■ 

Lemma 3.7 If gate §i satisfies that both of its inputs are outputs of other gates, 
then the effective contribution of JJ{ to the dimension of the range of S is at most 
9. 

Proof. The argument is the same as that for Lemma 3.6, except that the 
group acting independently on each input bit is used. This group has dimension 
7, so there are still 9 dimensions to parametrize the remainder of U gi in local 
neighborhood patches. ■ 

The dimension of the modified domain of f$ can now be counted as 
follows: All gates contribute 9 dimensions for a total of 96 dimensions. At most 
n inputs do not come from a previous gate. Those n inputs effectively contribute 
3ra + 1 dimensions, where the +1 comes from the phase of the first gate. 

The arguments given so far show that with b gates, unless 96 + 3n+ 1 > 
d, where d is the dimension of the range manifold, the measure of the range of 
fs, being a finite union of submanifolds of dimension at most 96 + 3ra + 1, is 0. 
This implies that 6(2, n; k) > |(2 n+1 - k) - \n - §. U 

4 Approximation of unitary operations by quantum circuits 

We begin by considering approximability by quantum circuits of unitary oper- 
ators in the || • 1 1 2 , A; norm. Let £7( n ,fc, 3 ,&, e ) be the set of unitary operators U such 
that there exists a circuit in Gg'' n ^ which computes a unitary operator V with 
\\U - V\\ 2 ,k < e- Let Exp(2,a;) = T and p(x) = x^l - x 2 /4. 

Theorem 4.1 For a > and (1 + a)e < a/2, 
»k(Un,k,g,b,e) < Ex P ( 2 , 2 4fl 6(log(6) + log(n) + log(2/(ae))) 
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k (T log(lM(l + a)e) - log(l/(l - (1 + a)V/2)))) 



The constants in the bound of Theorem 4.1 can be improved by more 
careful estimation. The main consequence of the theorem is the following corol- 
lary. 

Corollary 4.2 For fixed g and e < Ifb\og(b) = o(k2 n ) then fJ.k(U n ,k,g,b,e) = 
o(l). 

Note that y/2 is the expected distance of a random state to a fixed 
state, so that if 61og(6) = o(k2 n ), then for most \u) the circuits with less than b 
gates can do essentially no better at mapping |0) to \u) than a random unitary 
operator. 

Proof of Theorem 4.1. The proof is based on a covering argument. We first 
replace the unitary operators implemented by a single gate by a finite set of 
operators which is sufficiently dense. The number of distinct circuits is now 
finite. With the right choice of parameters, the volume of unitary operators 
that can be approximated is then given by the sum of the volume of suitable 
balls around each of the finite set of operators which can be represented exactly. 

We first compute the volume of a <5-ball for the || 1 1 2 , A; norm relative to 

Lemma 4.3 Let U be a unitary operator and B(U,S,k) = {V G U n \ \\V — 
?7| |2,fc < <*>}• Then 

/ , \2" +1 -l-2i 

*-i (SJl - S 2 /4) 
, k (B(U,S,k)) < (l-* 2 /2)- fc l P V2 n + i_i_ 2i 



< (l-£ 2 /2)" fc (6y/l-6 2 /4) 



k2 n 



Proof. For a logical formula (f>, let [(f>] = 1 if (f> is true, and [(f>] = otherwise. 
Let \uq), . . ., be the first k columns of U. We have 



fi k (B(U, S, k)) = J du k (u)[u G MB(U, 6, k))} 

< /^(vjlve^M^^-l))] 

dfii(\u))[ \u) _L v and \\u) — \uk-i)\ < 6] 
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< 




2k+l 



where we have used Lemma 2.4. The first inequality now follows by induction. 
The second one follows by observing that k < 2 n , the dimension of the space. ■ 

Lemma 4.4 For S > 0, there is a subset U nt $ of U n with 



such that for any U G U n , there exists V G U nt $ with \\U — V\\z < <*>• 

Proof. For U G U n , each entry Uij has norm at most 1. We can construct a 
set A of operators (not necessarily unitary) by selecting each of the 2 n matrix 



U — i?'s entries are within p of 0. Hence \ \U — B\\2 < \\U — B\\p < p2 n . Choose 
p2 n = S/2. Construct U nt $ by selecting for each A G A a nearest unitary matrix 
(according to the \ \ \\f norm). This increases the maximum distance by at most 
a factor of 2, so that our choice of p is good. The cardinality of the resulting 



Choose S = ^ in Lemma 4.4. If U G then by Lemma 4.4 



and 2.3, there exists V G G g g such that \\V — C^||2,fc ^ (1 + a ) e - Here we 
used G 9t s to denote the set of operations achievable by a single gate whose 
unitary operator is in U 9t $. The volume of U nt k, g ,b,e is therefore at most the sum 
of the volumes of the (1 + a)e balls (relative to || • ||2,fc) around each member 
of G^ b i a \ The relative volume of these balls is bounded in Lemma 4.3. Since 



\G g 4 h ' n) <(M^g)0) h we s et 

»k(Un,k, g ,b,e) < (m(S,g)^ (l-(l + a) 2 e 2 /2)- k p((l + a)e) k2 " 

< Exp (2, 6(2^(log(6) + log(2/(ae))) + log(n)) 

- k (T log(lM(l + a)e) - log(l/(l - (1 + «) 2 e 2 /2)))) • 





set is less than the stated bound. 
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From the point of view of quantum computation, a bound on approx- 
imability of density matrices on a subset of the bits is critical. Let U nt i t k, g ,b,t be 
the set of unitary operators U such that there exists a circuit in G g b ' n ^ which 
computes a unitary operator V with vi t k(U,V) < e. 

Theorem 4.5 For a > and 2(1 + a)e < 1, 

vi,k(Un,l,k,g,b,e) < Exp (2, &(2 4fl (log(6) + log(4/(ae))) + log(n)) 
- ^(loga^l + ^e))). 

Asymptotically, Theorem 4.5 has the following consequence: 

Corollary 4.6 For fixed g and e < \, If blog(b) = o(k2 l ) then Hi^U n ,l,k,g,b,t = 
o(l). 

Corollary 4.6 is not tight in the same way that Corollary 4.2 is. The 
minimum average total variation distance to a distribution on N points is | > | . 
This gap can be improved by more careful estimates in Lemma 2.5. 

Proof of Theorem 4.5. The proof follows the same outline as the proof of 
Theorem 4.1. The analog of Lemma 4.3 is as follows: 

Lemma 4.7 Let U be a unitary operator and B(U, 8, /, k) = {V G U n \ vi t k(V, U) < 
8}. Then 

^ k (B(U,8,l,k))<(2e) k ^- 1 \ 

By Lemma 2.1, we can also use Lemma 4.4 for the norm, taking 
care to compensate for the factor of 2 in the inequality of Lemma 2.1. Choose 
8 = ae/2b in Lemma 4.4. The same argument used in Theorem 4.1 gives 

vi,kU n ,,, k , g , b , £ < (m(8,g)^ (2(l + a)ef 2 '- 1 ) 

< Exp (2, 6(2 4fl (log(6) + log(4/(ae))) + log(n)) 
-^(loga^l + aOe))). 
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5 Classical approximation problems 

There are two classical approximation problems to be considered. 

Definition. Let / : D C 2^ — ► {0, 1} be a boolean decision problem and 
g : D — ► 0, 1 a probabilistic function. Then g (probabilistically) decides f with 
advantage q if for all b £ _D 

Prob(</(&) = /(&))>±(l + ±). 

If g decides / with advantage g, then f(b) can be obtained correctly with high 
probability by 0(q 2 ) independent evaluations of g. 

To probabilistically decide / with advantage q by a quantum circuit 
(that is by a unitary operator U £ U n i with n 1 > n) means that for b £ D, 
a measurement of the first bit of U\b,0) yields f(b) with probability at least 




Definition. Let / : D C 2^ 2^ be a function and g : D ^ 2^ a proba- 
bilistic function. Then g guesses f with advantage q if for all b £ D, 

Frob(g(b) = /(&)) > - 

A checking oracle for / is a function h which on input 6, b' decides 
whether f(b) = b'. If g guesses / with advantage q and we have a checking 
oracle for /, then we can obtain the values of / with high probability by 0{q) 
independent evaluations of g, checking each evaluation using the oracle. 

To guess / with advantage q by a operator U £ U n i with n' > n means 
that J2b' {f(b)b'\U\b{))\ 2 > K Note that the sum is the square amplitude of the 
projection of the result onto the space determined by setting the first n bits to 
fib)- 

We now show that there are decision and guessing problems which are 
difficult for quantum computation. To avoid trivial cases, we assume that g > 2 
and b > 2. We can also assume bg > n' — n. Otherwise some of the additional 
input qubits do not appear as an input to a gate and therefore do not contribute 
to the computation. 

Theorem 5.1 The fraction of decision problems with domain D decided by 
elements of G^g' n ^ with advantage q is at most 

Exp(2, (log(6) + log(g) + log(n))62 4fl+1 - |Z>|). 
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Proof. There are 2' B ' decision problems for domain D. Consider a fixed U G 
U n i. How many decision problems can be decided probabilistically by V G U n 
with \\V — U\\2 < 8 with advantage ql Let / be one decision problem with 
domain D solved by a V G U n i with \\V — U\\2 < 8. For \\W — U\\2 < 8 and 
b G D we have 

^(^|6)-y|6)) < 2\\W\b)-V\b)\\ 

= 2\(W- U)\b) + (U -V)\b)\ 

< 2 \(W-U)\b)\ + 2\(U-V)\b)\ 

< 48. 

Thus if 8 < ^, then W either does not decide any decision problem with domain 
D with advantage q or it decides / with advantage q. This implies that the 
number of decision problems that can be decided by members of G^g' n ^is at 
most the number of - balls required to cover G^g' n \ Using Lemmas 2.3 and 4.4 
with 8 = j^, this number is at most 

< Exp(2, (1 + log(6) + log(g) + log(n'))62 4fl ) 

< Exp(2, (log(6) + log(g) + log(n))62 4fl+1 ), 

where we have used the assumptions rather loosly to eliminate n' in favor of n. 

■ 

Corollary 5.2 For a sequence of domains D n C 2^ n \ there are decision prob- 
lems that cannot be decided by quantum circuits with advantage q unless (log(6) + 
log( g ) + log(n))6 = fi(|D„|). 

Note that for interesting domains, D n is exponential in n. 
We can obtain a similar result for guessing functions. 

Theorem 5.3 The fraction of functions f : D C 2^ — ► 2^ guessed with 
advantage q by operators in G g b ' n ^ is at most 

Exp(2, (log(6) + log(g) + log(n))62 4fl+1 - (n - log(4 9 ))|D|). 
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Proof. The number of functions considered is 2 n ' B '. We obtain a bound on 
the number of functions guessed with advantage q by operators V within 8 in 
2-norm of a fixed operator U £ U n i. For a vector u, let Pb(u) be the square 
amplitude of the projection of u onto the space spanned by the states \b, V). We 
have ^6 Pb(u) = M 2 . 

Let \\V - U\\ 2 < S. Let be D, \u) = U\b) and \v) = V\b). Let 

S(U,b) = {c | P c (\u))>l/4q}, 
S(V,b) = {c | P c (\v))>l/q}. 

If V guesses / with advantage q, then f(b) £ S(V,b) for each b £ D. We have 
- \v)\ < S. We would like to estimate \S(V,b)\ S(U,b)\. 

\\u)-\v)\ 2 = £P c (|«)-h>» 

C 

> £ Pc(k)-|^)) 
ce5(y,6)\5(C/,6) 

> £ 1 

ces(v,&)\s([/,&) y 
= \S(v,b)\S(U,b)\± 



where we used 



Pc(\u) ~ \v)fl 2 > \P c (\v)fl 2 - P c (\u)fl 2 \ 

1 



> 



If we choose S 2 < i, then |S(T, 6) \ S(£/, 6)| < 1, so that S(T, 6) C We 
can estimate S(U,b) by observing that 

i = IHII 2 

c 

ces(u,b) 

^ 4(7 

c£5(C/,6) ^ 
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which implies that |5(?7, 6)| < 4q. With this 6, the set of functions / with 
domain D guessed by any V with \\V — U\\2 < -^j^ satisfies f(b) £ S(U,b) for 
each b £ D, so its cardinality is at most (4g)' B '. 

The number of functions / : D C 2^ 2^ that are guessed with 
advantage g by members of G g b ' n ^ can now be bounded by 

< (4qj D \4^bf i3 n' b 

< Exp(2, (2 + log(6) + i log(g) + log(n'))62 4fl + log(4g) |D | ) 

< Exp(2, (log(6) + log(g) + log(n))62 4fl+1 + log(4 9 )|D|). 



Corollary 5.4 For a sequence of domains D n C 2^, i/iere are functions f : 
_D n C 2^ — ► 2^ no£ guessed with advantage q by any member of G^g' n ^ unless 

(log(6) + log(g) + log(n))6 = fi((n - log(4 9 ))|Ai|). 

Again, the interesting cases are where \D n \ is exponential in n. 

A Appendix 

Proof of Lemma 2.4. Assume first that m is the dimension of the full space. 
We can give an explicit integral for fj,(B u<t f~l S): 

fi(B Uje nS) = / #sin(^) 2m - 2 ° 



o ^ vr/ r(m-l/2) 

for e < y/2. The integral is obtained by parametrizing the spheres {v \ \v — u\ = 
<*> < e ? M = 1} by the angle (in real Euclidean space) between the v and u. We 
can write 

sin((i) = (5 rS ~2m-2 

#sin(^) 2m - 2 = 



o Jo \A — x 2 

/•S v 2m-2 

< 



< 



o VT^P 



(2m- l)Vl - <5 2 ' 
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Substituting the upper bound in terms of e gives 

(, \ 2m — 1 

MU u ' £ J - (2m-l)V5FT(m-l/2) (l-e 2 /2) 
Using the Sterling approximation we can estimate for m > 2 

y/2n(m-l) 1 < r(m) < ^/4tt(to - 1) (^T^)" 1 • 

Thus 

Tim) r- — / m-1 Y" 1 ' 1 

, V ; , , < A/2(m-l)/e 

r(m-l/2) " V v ;/ \ m -3/2j 

/ / 1 \-( m ~ 1 ) 

= ,/2(m-l)/e 1-— 

V v ;/ V 2(m-l)y 

To obtain good upper bounds, use the inequalities 1/(1 — x) < e x /(l — x 2 ) and 
fl - a;) m > fl - mx). 



T(m) < 



\/2(ra - 1) 7-4 s 

V v ; 1 - \ (Mm - 1 



r(m-l/2) " v v '1 - l/(4(m- 1))' 
Therefore 



Mi(5 u , £ n5) < 



( e yr^74) 



2m-l 



< 



y/(2m- l)vr(l - l/(4(m- 1)))(1 - e 2 /2) 

/ / \ 2m — 1 



V2m - 1(1 - e 2 /2) 

for m > 3 and e < \/2- 

To complete the proof of the lemma, we show that in the general case, 
B u ,t H S C fl 5, where f is the real projection of u onto the subspace 

generated by S. We can assume that v is non-zero, as otherwise for e < \/2, 
B u<t fl S = 0. Let • denote the real inner product. Let x £ _B Uje fl 5. Then 
u-x = v- x < Note that u, x and v are in the same halfspaces determined 

by the hyperplane perpendicular to v and the one perpendicular to x (otherwise 
\u — x\ > -v/2)- Let cj) and tjj be the angles between u and x and between f and 
x, respectively. The inequality between the inner products means that (f> > if). 
Since u, x and v are unit vectors, \u — x\ > \v — x\, which gives the result. ■ 

Proof of Lemma 2.5. Let v £ A(N). Let H be the hyperplane which contains 
A(N). Let B be the e ball in the L\ metric around v. To compute ji(B fl H) we 



21 



can shift v to the origin and consider each orthant separately. After projection 
onto all but one of the negative coordinates, the intersection of B n H with an 
orthant with k positive coordinates is equivalent to the set 

{k N-l k "I 

x | x t > 0,^2 x t < e/2, ^2 x i<^2 x i\- 
i=l i=k+l i=l ) 

The total projected volume of the orthants is therfore given by 

KBnH) - \V)L ^^/jt* w^kj 

y( N ) I u /2 f-i 

^ [k ) (N - l)T(k)T(N - k) [ 1 ) 
1 ^ {N\ {N-2\, .^n^ 



< i (2e) N -\ 



Since the projected volume of A(iV) is p^y) , the result follows. ■ 

Lemma A.l Let X and Y be n X m operators such that X*X = Y*Y . Then 
there exists a unitary operator U such that UX = Y . 

Proof. Let Y = WDV be the singular value decomposition of Y, with D a 
non-negative diagonal matrix. By a change of coordinates, we can assume that 
V = I. This gives Y*Y = D*D. X*X = D*D implies that the columns of X 
are orthogonal, of lengths corresponding to the entries of D. Hence there is a 
unitary transformation W such that W'X = D. Letting U = WW 1 gives the 
lemma. ■ 
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